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1 Introduction 


Recently there has been spectacnlar progress on exact solntions for strongly interact¬ 
ing stochastic particle systems in one spatial dimension. Very ronghly, these are two- 
dimensional held theories with some similarity to two-dimensional systems of eqnilibrinm 
statistical mechanics and 1-1-1 dimensional qnantnm held theories, inclnding qnantnm 
spin chains. In these stochastic models the static correlations decay exponentially, bnt, 
becanse of conservation laws, the dynamic behavior shows non-trivial scaling properties. 
The physical interest in snch models resnlts from being representatives of hnge nniver- 
sality classes. This is in direct analogy to the two-dimensional Ising model. Its exact 
sointion is restricted to nearest neighbor conplings and zero external magnetic held. Bnt 
the thereby predicted critical behavior is representative for a mnch larger class of models, 
inclnding real two-dimensional ferromagnets to which, at hrst sight, the Ising model wonld 
be a very crnde approximation. 

Onr review attempts to answer “What is cnrrently known abont the one-dimensional 
Kadar-Parisi-Zhang eqnation?”. This will cover mathematically rigorons resnlts, bnt also 
important conjectnres and intrigning replica compntations. In addition, we discuss dis¬ 
crete models which in a particular scaling limit converge to the KPZ equation. They 
constitute the universality class of the KPZ equation alluded to in the title. Such models 
have been, and still are, an important tool to arrive at information on the KPZ equation 
itself. Our program has obvious disadvantages. It mostly ignores in which temporal order 
discoveries were made. It also misses largely the amazing progress on other models in 
the KPZ universality class. Fortunately there are excellent reviews available [1-8], which 
cover further aspects of the story. 

1.1 Interface motion. In our article we focus on the KPZ universality class for growing 
fronts and interfaces in two-dimensional bulk systems. With the hindsight of theoreti¬ 
cal activities over many years, let us hrst recall the physical set-up under consideration. 
One considers a thin him, for which two bulk phases can be realized such that they are 
in contact along an interfacial curve. One of the phases is stable, while the other one 
is metastable with a very long life time. The dynamics of the bulk phases has no con¬ 
servation laws and is rapidly space-time mixing. Therefore, while huctuating, the bulk 
phases are not changing in time and the only observable motion is at the interface, where 
the transition from metastable to stable is fast, while the reverse transition is very much 
suppressed. Thereby the stable phase is expanding into the metastable one, accompanied 
by characteristic interfacial fluctuations. If both phases would be stable, then there is no 
net motion, and the interface fluctuations are of size with Gaussian statistics on a 
sufficiently coarse scale. In the KPZ set-up on the other hand, in the frame where the bulk 
phases are at rest, the interface has a non-zero net velocity and fluctuations are of size 
a phenomenon known as kinetic roughening. There is no hne tuning of parameters. 
Any asymmetry in the growth process will lead to the scaling. Such motion comes 
with interfacial fluctuations which have non-Gaussian statistics, hence the prime interest 
in exact solutions. 

The most prominent physical realization is a thin him of turbulent liquid crystal [9,10]. 
In fact, the bulk phases are driven by a time oscillating, spatially uniform electric held. 
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Figure 1; The bold line is the current height. The transition from V to A occurs with 
rate q, from A to V with rate p, q > p. 


and thus not in thermal equilibrium. At a carefully chosen point in the phase diagram, the 
DSM2 phase (stable) is coexisting with the DSMl phase (metastable). In the experiment, 
a point seed of the stable phase is planted in the bulk metastable phase and the growth 
of the DSM2 cluster is observed. In further series of experiments also the case of a stable 
line seed is realized. Another example is Eden cluster growth. Here one starts from a 
point seed (stable phase) and randomly adds extra material to the surface of the already 
existing cluster. Hence the metastable phase is just “all empty” and the precise ballistic 
deposition mechanism is not specihed by the model. On a large scale, the cluster has an 
approximately circular shape. But there are substantial shape fluctuations, which have 
the same statistics as observed for the turbulent liquid crystal [11,12]. 

Is there a model system which captures the essential features and is still amenable to 
a theoretical analysis? The arguable simplest model is the two-dimensional ferromagnetic 
Ising model with Glauber dynamics at zero temperature and a suitably scaled positive 
magnetic held. The -|- phase is stable and the — phase is metastable. We assume that the 
-|- phase is bordered from above by a lattice height function, h{j, t), which is a broken line 
with slope ±1, see Fig. 1. The admissible growth sites are the local minima of h, at which 
independently a — spin hips to -|- spin after an exponentially distributed waiting time, the 
reverse process occurring with a much smaller rate. In principle, clusters of the “wrong” 
sign could be spontaneously generated inside the bulk phases. But at low temperatures 
this is a very unlikely event and is completely suppressed at zero temperature. If instead 
of the height one considers the local slope with value ±1, then the interface dynamics 
amounts to a random exchange of a neighboring —h pair, the reverse process being less 
likely. This Markov jump process does not satisfy detailed balance, which is one reason 
for the unexpected dynamical behavior. The one-sided growth limit, g = 1, is known as 
TASEP, which has been a useful tool in obtaining exact solutions, see [3]. 

Kardar, Parisi, and Zhang [13] pose the same question, but proceed along a very 
different route, motivated by the Ginzburg-Landau theory of critical dynamics. Firstly, 
on a mesoscopic scale the interface is modeled as the graph of a height function h{x, t) 
with time t > 0 and space x G M. In the case of contact between two stable phases one 
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writes down the Langevin equation 

dth = udlh + . ( 1 . 1 ) 

Here u is the surface tension and \/D^ is a noise term of strength \/D describing the 
random back and forth between the two phases. ^ {x, t) is modeled as space-time white 
noise, which is Gaussian with mean zero and covariance 

{^{x,t)^{x',t')) = 6{x - x')6{t -t'). (1.2) 

On the basis of (1.1) one hnds that, if initially flat, the interface develops fluctuations of 
size 

On the other hand, if the interface is stable-metastable, the net motion is modeled by 
adding to (1.1) a nonlinear drift, which is allowed to depend only on d^h so to maintain 
the invariance under the global shift of h to h -|- a. For isotropic growth in the plane the 
nonlinearity can be easily guessed. If locally there is a linear piece, h{x, t) = a + 'dx, 
then h{x,t + dt) — h{x,t) = noo(l + = F(id)dt with Voo the asymptotic radial 

growth velocity. The precise form of the nonlinearity should not matter on large scales 
and, expanding at = 0, one arrives at the one-dimensional KPZ equation 

dth = lX{dxhY + vdlh + (1.3) 

with A = F"(0). In principle, we should keep the strength of nonlinearity, diffusion, and 
white noise as free parameters, but by rescaling x, f, h the conventional choice A = H = 1, 
z/ = or any other, can always be achieved. Also the sign of the nonlinearity plays no 
role as can be seen from the switch of h to —h. 

1.2 Cole-Hopf solution. The continuum description (1.3) of the interface has its math¬ 
ematical price. The solution to the linear equation (1.1) is locally a Brownian motion in 
X, i.e. it has everywhere cusp singularities as ±^/\x\. Hence dxh of (1.3) is a distribution 
whose square does not easily make sense. Fortunately there is a simple and instructive 
way out, which however is special for the quadratic nonlinearity. To keep things simple, 
set X = 1 = D, u = 1/2 and introduce the Cole-Hopf transformation 

Z{x,t) = (1.4) 


Then Z is the solution of the stochastic heat equation, 

dtZix^t) — \dlz{x, t) + t). 


(1.5) 


The additive noise in (1.3) has been turned into multiplicative noise. Using the Feynman- 
Kac formula, the “solution” to (1.5) with initial condition Z{x,0) can be written as 


Z(x,fj = [e-fi''* «*<*>'*>Z(6(0),0) 


( 1 . 6 ) 


Here refers to the expectation w.r.t. an auxiliary standard Brownian motion b{s), 0 < 
s < t, with b(t) = X. Despite b{s) being continuous, as written, the integral in the 
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exponential is ill-defined, since is too rough. Still we expand the exponential and 

carry out the Brownian motion average with the result 


OO n n 

Z{x,t) = ^ dy 

n=0 do 




I0<ti<...<t„<t 


n « 

^n(/ dxj^{xj,tj)^pt{y,0;xi,ti,...,Xn,tn;x,t)Z{y,0). (1.7) 


i=i 


Here pt is the time-ordered multi-time transition probability of Brownian motion. Eq. 
(1.7) is a multiple stochastic integral, and we choose the forward (ltd) discretization in 
every tj. Then each summand of (1.7) is a well dehned random variable, the sum converges 
in T^(^), and Z{x,t) > 0. We thus define the (Cole-Hopf) solution to the KPZ equation 
as 

h(x,t) = \ogZ{x,t). (1.8) 

While the sum (1.7) seems to be useful only for studying the short time behavior of the 
solution, it does provide us with a well-defined starting point. 

The singular structure of the KPZ equation comes from insisting on spatial white noise. 
If one would replace ^(x, t) by the smoothed version ^^{x, t) = J (p{x — x')^{x', t)dx' with 
(p > 0, even, smooth, of rapid decay at inhnity, and normalized to 1, then the solution to 
(1.3) is well-dehned and given by (1.8) with the random partition function as dehned in 
(1.6). We now introduce an ultraviolet cut-off of size e by choosing the ^-function sequence 
9 ?e(x) = e~^Lp{e~^x) and substituting the noise ^ by its smoothed version Such 

solution is denoted by hfix, t). As established some time ago [14,15], in the limit e —)■ 0 one 
recovers the Cole-Hopf solution, provided one switches to a moving frame as hfix, t) — vfi 
with velocity = e~^ J (p{xfidx diverging in the limit. In held theory language this 
amounts to an inhnite energy renormalization as familiar from P{ 4>)2 quantum held theory. 

In this context Hairer recently developed a solution theory, for which he was awarded 
the 2014 Fields Medal. His theory works for the KPZ equation, as well as a large class 
of other singular stochastic partial diherential equations, and for general cutohs [16,17]. 
For a specihc cut-oh, which does not pass nicely through the Cole-Hopf transformation, 
a somewhat more straightforward analysis can be found in [18]. Of course, the solution 
theory studies the small scale structure of solutions, and not the large scale, at which the 
universal behavior of interest is observed. 


1.3 Directed polymers in a random medium. Before diving into further details, let 
us recall an important observation made already in the 1986 paper. To be specihc, let us 
choose the initial condition Z{x, 0) = 6{x). Then b{s), 0 < s < f in (1.6) can be regarded 
as a directed polymer starting at time 0 at 0 and ending at time t at x. It has been 
assigned the random energy 

ds^{b{s),s) (1.9) 

as the sum over the potential ^ along the path b{s). With this perspective, Z{x,t) is a 
partition function, the sum being over all directed polymers. Z{x,t) is random, because 
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of the random potential Thns we really stndy a disordered system from equilibrinm 
statistical mechanics. The height of the interface, the object of physical interest, is simply 
the (random) free energy. To leading order it is proportional to t, which corresponds to 
the non-random displacement of the height. Thus physically our goal is to determine the 
fluctuations of the free energy. In the theory of spin glasses such hne details are rarely 
studied. 

1.4 Invariant measures and scaling. If one introduces the slope u = dxh, then u 
satishes the conservation law 


dtu + dx{ — — udxU — = 0 . 


( 1 . 10 ) 


The inviscid Burgers equation corresponds to the case u = 0, D = 0. Eq. (1.10) is 
hence referred to as noisy, or stochastic. Burgers equation. The linear equation, A = 
0, dehnes a Gaussian process. It has a one-parameter family of stationary measures, 
which are labeled by the mean {u{x)) = 'd and otherwise are Gaussian white noise with 
covariance {u{x)u{x')) — = a'^6{x — x'), a = \/ D/2v. Surprisingly enough, adding the 

nonlinearity does not modify the invariant measures. For this property to hold one only 
has to check that Gaussian white noise is invariant under the flow generated by solution to 
the deterministic evolution equation dtU = ^^XdxV?- Formally its vector held is divergence 
free and the invariance follows from considering the time change of the action. 


d 

dt 


u^dx = A 


udxU^dx = |A 


drU^dx = 0 . 


( 1 . 11 ) 


Our result can be translated to the KPZ equation [19]. One chooses as initial conditions 
the two-sided Brownian motion aB{x) + 'dx with arbitrary slope ■d, i.e. 5(0) = 0 and 
{B{x),x > 0}, {B{x),x < 0} are two independent standard Brownian motions. Then the 
solution h{x, t) at hxed time t is again two-sided Brownian motion with slope 'd except for 
a random displacement, h{x,t) — h{0,t) = aB{x) X-dx in distribution with B{x) another 
two-sided Brownian motion. Note that h{0,t) and B{x) are not independent. 

With this observation one easily determines the scale relevant for huctuations. We 
consider the solution h{x,t) to (1.3) and dehne the rescaled huctuation held 

h^{x, t) = e^h{e~^x, e~^t ). (1-12) 


Using the scale invariance of white noise yields 

dtK = \\e^-"-\dxKf + (1-13) 

To keep the stationary measure unchanged forces one to have the huctuation exponent 
b = 1/2. Thus to maintain the nonlinearity of order one requires 

z = 3/2. (1.14) 

Our rescaled equation takes the form 

dth = lX{dxhf + • (1-15) 
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This scaling h e t ~ e x ~ e ^ is more commonly thought of as h rsj 
X ~ This is the scale in which one observes the non-trivial fluctuation behavior 

throughout the KPZ universality class. As to be noted, to extract from the KPZ equation 
universal, model independent properties one still has to study its long time limit. 

To be more precise than their size, one must take into account the initial conditions. 
Physically natural options are singled out. One is droplet growth, e.g. the Eden cluster, 
for which the deterministic shape has a local curvature decaying as 1/t. Another natural 
choice would be growth on a flat substrate, modeled as that h{x, 0) = 0. As a consequence 
h{x,t) will stay flat on average. One also could imagine to evolve for a long time span r 
and then consider height differences as h{x, t + t) — h{x, r), which means that the initial 
height statistics has become stationary and one would have to compare with the KPZ 
equation for stationary initial conditions. 

As will be discussed in more detail in Sect. 2, for the curved case there is an exact 
solution for the distribution of h{0,t) for alH > 0 [20,21]. From it one deduces that 

h{0,t) ~ Voot + GUE (1-16) 

for large f, where v^o is the asymptotic growth velocity and P = ||A|(T'^. The random 
amplitude Cgue has the distribution function 

IP(Cgue < s) = det(l - PsKaiPs) ■ (1-17) 

The right hand side involves a Fredholm determinant over (M, dx) , see [22], Section 17, 
for Fredholm determinants. Pg is the projection operator onto the interval [s, oo) and Kai 
the Airy kernel 

poo 

^Ai{x,y) = / duAi{x + u)Ai{y + u) (1-18) 

Jo 

with Ai the standard Airy function. For the transverse exponent 2/3, sharp bounds are 
reported in [23]. For stationary initial conditions an exact expression for the distribution 
of h{0,t) has been recently accomplished [24] and an asymptotics corresponding to (1.16) 
can be deduced. 

For the sharp wedge also a replica solution has been achieved [25]. By the same method 
stationary initial conditions can be handled [26,27], obtained prior to [24]. However for 
flat initial conditions currently we have to rely on a replica solution only [28,29]. Its 
structure is much more complicated than in the other two cases. The replica solution 
tends to be supported by mathematically rigorous computations on discrete lattice type 
models in the KPZ class [30,31]. What is more certain from such discrete models is that 
for large t one should have 

/l(t,0) ~ Uoo^ + (rf)^/^CGOE , (1-19) 

where the random amplitude Cgoe has the distribution function 

IP(Cgoe < 2s) = det(l — PqBsPq) (1-20) 

with kernel Bs{x,y) = Ai(x + ?/ + s), see [32,33] for this particular representation. 
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The distribution functions (1.17), (1-20) were obtained first by Tracy and Widom [34] 
in apparently completely unrelated problems, namely as the distribution of the largest 
eigenvalues of GUE/GOE random matrices. Just to briefly recall the GUE case, one con¬ 
siders NxN complex Hermitian matrices. A, with the distribution Z~^ exp[—(2A^)“^try4^](iy4. 
A has the real eigenvalues Ai < ... < Xn, where the factor N~^ in front of tiA'^ ensures 
that their typical distance is of order 1. For the largest eigenvalue it is proved that 

Xn^2N + N^^^Cgue ( 1 . 21 ) 

for large N. Totally unexpected, the same random amplitude as for sharp wedge KPZ 
growth appears, compare with Eq. (1.16). The GOE case is in complete analogy with real 
Gaussian symmetric matrices. The GOE/GUE Tracy-Widom probability distributions 
can also be represented in terms of special solutions of the Painleve II equation, indicating 
a link to classical integrable systems. 

To return to the numerical/experimental side, if one wants to compare the KPZ pre¬ 
dictions (1.16), (1.19) with the data, one hrst has to carefully determine the value of P. 

A is measured through the second derivative of the slope-dependent growth velocity. For 
example, in the case of isotropic growth, one would have F{dxh) = noo\/l + {dxhy and 
hence A = v^o- The coefficient is more difficult to extract and involves larger uncertain¬ 
ties. In [9] it is used that = lima;^oo \x\~^Va.T{h{x,t) — The sampled 

distribution for the height at origin at various times plotted in these units should fall on 
top of each other and agree with the respective Tracy-Widom distributions. 

1.5 Outline. We first discuss the exact solution with sharp wedge initial data and then 
address in more detail the link of lattice type models for interface dynamics to the KPZ 
equation. The second main item is the replica solution. While this method is forced 
to work with diverging sums, it is still of interest as a tool and in connecting to exten¬ 
sive work on spin glasses and quantum quenches. We will also explain a closely related 
approach which deals only with convergent sums. 


2 Exact solution of the KPZ equation for initial sharp 
wedge 


On the KPZ level, droplet growth is induced by the sharp wedge initial condition 


h{x, 0) = —(5 ^ |a;| — log(2J), J —)■ 0 . 


( 2 . 1 ) 


The solution develops immediately a parabolic shape with superimposed fluctuations, 
which thus is taken to represent the top part of a droplet. The corresponding initial 
condition for the stochastic heat equation reads Z(x,0) = 6{x). We note that 


Z{x, t) = 


Z{x, t) 


{Z{x,t)) 

equals the partition function relative to the Brownian bridge 


w{s) = b(s) - j(b(t) - x) , 


( 2 . 2 ) 


( 2 . 3 ) 
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Figure 2: Probability density function -F^(s) for time t from short times (red, lower curves) 
to long times (green, upper curves) for t ranging from 0.25 to 20,000, see [36] for further 
details. For t ^ oo, the density converges to the GUE Tracy-Widom distribution Fque 
( upper blue curve) and, for f —)■ 0, F)'(s) becomes a Gaussian (rightmost blue curve) with 
mean and variance increasing respectively as and on the scale set by the 

Tracy-Widom distribution. 


0 < s < t. Expectation w.r.t. the Brownian bridge is denoted by Ebb- Then 


Z(x-|-a, f) = Ebb ( exp J ds^{w{s) + aj, s) 


(2.4) 


Since ^ is white noise, Z{x + a, t) has the same distribution as Z{x, t) which implies that 
X e-)■ Z{x^t) is a stationary process. On the level of the height function this means that 
h(x, t) +x‘^/2t is stationary in x. In fact the exact solution has still a shift by —Thus 
the object of interest is the deviation, rj, from the deterministic downwards parabolic 
prohle as 

r]{x,t) = h{x,t) + j^x"^ + , (2.5) 

where h{x, t) is the Gole-Hopf solution of the KPZ equation with sharp wedge initial data. 
By stationarity it suffices to consider the distribution of ri{0,t). 

The exact solution was obtained a few years ago [20,21] and provides an expression 
for the generating function. 


( exp [ - exp(? 7 ( 0 , t) - yjs)] ) = det(l - PoKt,sPo ), (2.6) 

valid for all real s G M and f > 0, where 

7t = (t/2)'/^ 
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The right hand side involves a Fredholm determinant over dx), Pq projects on M+, 

and Kt,s is the smoothed Airy kernel 

^ + “)Ai(l/ + u). (2.7) 

Note that the generating function (2.6) determines uniquely the distribution function 
Pf(s) = P(r7(0,t) < s). There is an explicit inversion formula yielding Pt(s), again involv¬ 
ing Fredholm determinants [20]. Numerically such determinants can be computed very 
efficiently through approximating the integral kernel of Kf s by a hnite-dimensional matrix 
of the form Kt,s{xi,Xj)\ij=i^,„^n with carefully chosen base points {xj,j = [35]. 

In our application n = 40 suffices already to achieve a good precision. In Fig. 2 we 
display the probability densities P/(s) in dependence on t. For short times the required 
matrix dimension n increases and the numerics becomes slow. In that limit one returns 
to the sum (1.7) to obtain that the initial fluctuations are Gaussian with a width of order 
At later times a characteristic peak develops which overshoots somewhat to the 
left and hnally settles, on the scale at the GUE Tracy-Widom distribution. The 
formula for the limit distribution is deduced from (2.6) by substituting ri{x,t) — by 
~ ^)- Then, taking the limit f —)■ cxo on both sides, one obtains 

lim P( 7 j"^r 7 ( 0 , t) < s) = det(l - PsKaiPs) ■ (2.8) 


From the exact solution one deduces also some information on hnite time corrections. 
The shape of the distribution is fairly rapidly attained with the slowest mode being the 
mean which approaches its limit as Similar behavior is reported from numerical 

simulations and from the experiment, with a small twist. As can be seen from Fig. 2, for 
the KPZ equation the hnite time correction has a negative sign, while in the other cases 
the sign is positive. This teaches us that, while the power is apparently universal, 
the prefactor depends on the particular model [37]. 

The obvious next step would be to consider the joint distribution of y(0,f),y(x,f). 
No closed expression as in (2.6) is available. But, based on universality, the long time 
behavior is conjectured with a high level of conhdence. For simplicity let us consider the 
joint distribution of 7j~^'^^?7(0, f), f). We know that these random variables stay 

of order 1 as t —)■ oo. For hxed x, they will be perfectly correlated as f —)■ oo, while 
they will become independent if simultaneously |a:| —)■ oo sufhciently fast. As anticipated 
already in the line below (1.15), the interesting scale is x = 

Conjecture. For every w it holds that 

lim P(7"^/^r7(0, t) < Si, t) < Sa) = det(l - ), (2.9) 


/oj 

where Ksi\s 2 ;w is the extended Airy kernel. 

The determinant needs some explanations. acts on L^(M x {1,2}) and is given by 


/^( 2 ) 

si,; 


fPs, 0 W Km -e-l-l^(l - Km) 
V 0 PsJ \e\^\^KM Km 


( 2 . 10 ) 
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Ps is the projection onto [s, oo) and H is the Airy operator, 

( 2 . 11 ) 

Then Kx\ is the projection onto all negative eigenstates of H, K^i = P{{H < 0}), hence 
e-l'^l'^(l — Kai) and are bounded operators. More generally one expects that 

tc 1-4 ^^^ri{2'y^w,t) as a stochastic process converges to the Airy process [38], dehned 

through hnite dimensional distributions given by expressions generalizing (2.9) to multi¬ 
space points. 

With Eq. (2.9) one has a tool to study the asymptotic height-height correlation 

lim ((7r^^^^(0,t) - = g{w). (2.12) 

t^oo 

g{w) ~ |tc| for small |tc| and 

g(w) ~ — 2w~‘^ (2.13) 

for |t<;| —)■ cxo with the variance of the GUE Tracy-Widom distribution, see [39]. For 
hner properties one has to rely on the numerical evaluation of the Fredholm determinant 
in (2.10). 

Correlations at the same location but two different times, as e.g. the joint distribution 
of have been regarded as inaccessible by current methods. But, 

based on the zero temperature semi-discrete directed polymer, an asymptotic expression 
has just been posted [40]. Prior it has been argued that under minimal assumptions this 
joint distribution remains non-degenerate in the limit t ^ oo [41]. 

The proof of (2.6) is fairly indirect. The general strategy is to start from a discrete 
growth model, for which sufficient details on the solution with sharp wedge initial data 
are available. One takes a suitable scaling limit in which the discrete model is known to 
converge to the KPZ equation. Then the limit of the generating function for the discrete 
growth model will yield an identity similar to (2.6). For the sharp wedge one convenient 
approximating model is the ASEP with step initial conditions [15], see Sect. 3.2 for 
more details. Tracy and Widom [42] derived a Fredholm determinant for the distribution 
of the m-th particle at time t, which is sufficiently concise for an asymptotic analysis. 
A similar strategy has been implemented very recently for stationary initial conditions, 
h{x, 0) = B{x) with B{x) two-sided Brownian motion [24]. In this case the approximating 
model is constructed through the theory of Macdonald processes. 

For flat initial conditions, h{x, 0) = 0, the situation is mixed. For PNG, TASEP, 
and the zero temperature semi-discrete directed polymer the large time statistics starting 
from flat initial conditions have been investigated [43-46]. Since these models do not 
have an asymmetry parameter, no information on the KPZ equation can be extracted. 
Currently the only candidate is the asymmetric simple exclusion process [30,31]. However 
the naturally arising generating functions are not so well suited for an asymptotic analysis. 
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3 Approximating the KPZ equation 

For models close to a realistic physical description there is always a smallest length scale, 
say the size of the smallest object attached at, resp. removed from, the interface. Thns 
the issne arises how snch models are related to the KPZ eqnation. In textbook discussions 
the focus is on the expansion 

+ {dxh)‘^ = 1 + \{dxhy‘ + 0{{dxhY ), (3-1) 

which looks like a small slope expansion. But h{x) has locally the statistics of a Brownian 
motion and d^h has many spikes. So perhaps one should hrst average over some inter¬ 
mediate scale. Through a deeper study of a few model systems, a perspective different 
from Taylor expansion has been developed. Namely the KPZ equation is derived from a 
microscopic model in the limit of weak drive. We illustrate with two examples, but hrst 
take up the issue of nonlinearities different from quadratic. 

3.1 Viscous random Hamilton-Jacobi equations. These are stochastic partial dif¬ 
ferential equations of the form 

dth = F{dxh) + vd^h + (3.2) 

with a general non-linearity F. is space-time white noise smoothed at scale e, see the 
paragraph below (1.8), and the problem is to make sense of the solution as e tends to 
zero. A natural example of such a non-linearity would be the lateral growth mechanism 
with F('d) = \/l -|- as in (3.1). In order to approximate (3.2) by the KPZ equation, 
what would be done is to expand 

F(d) = F(0) + F'(0)d + iF"(0)d2 + • • • . (3.3) 

The K(0) and K'(0) can be trivially removed by height and spatial shifts. Arguing that 
the higher than quadratic terms are negligible one obtains in the new frame the KPZ 
equation (1.3) with A = F"(0). 

A more convincing approach is to take the ultraviolet cut-off e —>■ 0 simultaneously 
with a weak nonlinearity of size y/e. Rescaling as in (1.12) with h = 1/2^ z = 2 leads to 

dtK = e~^F{e^/‘^dxh^) + (3.4) 

Recall that the stationary measures of the KPZ equation are of the form 

h{x) = aB{x) + 'dx , (3.5) 

where B{x) is a two-sided Brownian motion. The stationary measures of (3.4) are ex¬ 
pected to be approximately of the same form, but with B smoothed out on scale e. One 
expects that h^{x,t) viewed near x is close to a stationary measure with the appropriate 
slope. The rapid fluctuations replace F('d) by .F('d), F denoting average with respect to 
such a locally stationary measure. Then, for small e, the effective equation reads 

dthe = e~^F{e^/‘^dJk^) + i^dlh^ + . (3.6) 
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is less wiggly than and it is now allowed to take the formal limit to arrive at the KPZ 
equation (1.3) with A = .F"(0) instead of F"(0). Such a claim is proved for polynomial 
nonlinearities in [47]. 

As an example, one may want to make sense of the quartic KPZ equation 

dth = {dxhY + T^dlh + ^ (3.7) 


by smoothing out the noise and rescaling so that all terms are of unit order in the limit. 
Then F('d) = ((uXc + "d)^) with Xq a standard Gaussian and one merely arrives back 
at the quadratic KPZ equation with a renormalized A = .F"(0) = 12 (T^(Xq) ^ 0. We 
learn that the quadratic KPZ equation is the only case in which it makes sense to have 
delta-correlated noise. 

3.2 Asymmetric simple exclusion process. The asymmetric simple exclusion pro¬ 
cess models particles on a one-dimensional lattice Z performing nearest neighbor random 
walks, with an asymmetric jump rate and subject to the exclusion rule of at most one 
particle per site. The jump rates are p G [0,1] for jumps to the right, i.e. from j to j -|- 1, 
and q = 1—p for jumps to the left, all jump trials are independent, and jumps attempted 
onto an already occupied site are suppressed. The ASEP height function is 

defined by assigning slope —1 to an empty site and slope 1 to an occupied site, compare 
with Fig. 1. t) evolves with its own Markovian dynamics, which flips V ^ A at 

rate q and A -w v at rate p. For q > p the general trend is an upwards moving height 
function. 

In this model there is a remarkable discrete version of the Cole-Hopf transformation 
(1.4) as first noted in [48]. We define the partition function 


_ g-(p+g-2y^)ig|(logT)hASBPQ^j) ^ 


T = p/q. 


(3.8) 


It satisfies the stochastic differential-difference equation 

^^ASEP _ . (3.9) 

Here A/(x) = f{x -l- 1) — 2/(x) -|- f{x — 1) is the discrete Laplacian and dM is a family 
of derivative martingales with {dMj, dMji) = 0 for j ^ j'. As in the previous example we 
adopt the weakly asymmetric limit 


q — p = . (3.10) 

Then \\ogT ~ at diffusive scales, time = e~H, j = e~^x, Eq. (3.9) 

converges to the stochastic heat equation (1.5). Taking logarithms one has established 
that the motion of the discrete ASEP height function converges to the Cole-Hopf solution 
of the KPZ equation [15]. 

Let us specialize from general initial data to step initial conditions, i.e. 0) = 

— |j|. Then to leading order h^®^^(0,t) ~ —For one-sided jumps, at times so short 
that the discrete nature of h^®^^(0,f) is still resolved, numerically one observes a discrete 
probability measure sitting on top of the Tracy-Widom probability density except for a 
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global shift of order The modeling through the KPZ equation is not so meaningful. 

Only at sufficiently long times the distributions from both models coincide. On the other 
hand at sufficiently weak drive there is an intermediate scale on which ASEP and KPZ 
agree very well. At this point one should worry whether in the limit of weak drive the 
universal asymptotics might be lost. A famous example are the long time tails of current 
correlations in fluids which decay as in the physical dimension. On the other hand 
kinetic theory, a low density approximation, predicts exponential decay. In our context, 
the exact solution conhrms that the KPZ equation just barely maintains the nonlinear 
effects. 

3.3 Directed polymer in a random medium. A further instructive example is the 
directed polymer representation as in (1.6). Probabilistically it is rather obvious how 
to discretize space-time as 1?. b{t) is turned into a discrete time random walk Sn with 
n G N, Ao = 0, and independent increments — S'n = ±1 with probability White 
noise is replaced by the family of i.i.d. random variables G Z. Let Sn be an N 


step random walk. We assign an “energy” by 

N 

(3.11) 

n=l 

Then the point-to-point partition function reads 

Z°^P(iV,M)= 2-^exp[pE^{S)]. (3.12) 

We require that < oo at least for small jS. M) satishes a one-step local 

recursion relation, 

+ 1, M) = M + 1) + M - 1)) . (3.13) 

Taking logarithms, with h{n,M) = log Z^^^{n, M), the recursion becomes 

h{n + 1, M) = log + 1, M), (3.14) 

which can be viewed as a particular growth process with additive i.i.d. noise. By analogy 
one expects, that 

logZ°^P(iV,0) ~ v^N + iTN)^/^CGVE (3.15) 


by KPZ universality. \min^oo'n-~^logZ^^^{n,Q) = Uoo a.s. is a consequence of subad¬ 
ditivity. The asymptotic velocity Voo can be computed only for few models [49]. Our 
understanding of the fluctuations is restricted to the case of C being log-Gamma dis¬ 
tributed [50]. 

To derive the Cole-Hopf representation we choose a lattice spacing e. Brownian motion 
limit then requires n = [-J denoting integer part. To have the energy (3.11) to 

converge to the continuum energy from (1.6) one has to scale the noise weak as ^/e. This 
leads to the intermediate disorder scaling, 

e"""'^(("/^)'/')%/2)i/2(Le"2fJ, le-^x \), (3.16) 
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where A(/9) = E(e^^). This partition function converges to the solution of the stochastic 
heat equation with Z{x, 0) = 5{x) [51]. The factor 2“^/^ simply comes from the fact that 
(3.13) lives on the even sites of the two dimensional lattice. It should be emphasized 
that the result is independent of the distribution of the randomness, only relying on 6 
hnite moments of min(0,(^). Again, taking logarithms we see that in these models the 
free energy converges to the respective solution of the KPZ equation in the intermediate 
disorder limit. 


4 Bethe ansatz and replica solution 

As noted already early on [52], there are closed evolution equations for the moments 
of Z{x,t) as dehned in (1.6) . Let us consider the A^-th moment, now written with N 
auxiliary Brownian motions, called the replicas. Denoting the white noise average by (•), 
one arrives at 

N N 

(nzfe.t))=(■‘•1) 

i=i i=i 

The white noise average can be carried out explicitly and is given by the exponential of 



In this form, the Ito discretization amounts to omitting the diagonal term i = j. The 
double time integration reduces trivially to a single one. Hence, using Feynman-Kac 
backwards, one obtains 


N 

{l[Z{xj,t)) = {xi,...,XN\e-^^'\Z^^). (4.3) 

1=1 

Here is the A^-particle Lieb-Liniger quantum Hamiltonian on the real line with at¬ 
tractive (5-interaction, 


N N 

1=1 *^ 1=1 

The quantum propagator acts on the product wave function denoted by 

and is evaluated at the point (xi, ...,xn)- Since Z^^ is symmetric, only the restric¬ 
tion of exp [—LfAft] to the bosonic subspace in L^(M^) is required. As a result the right 
hand side of (4.3) is a symmetric function, as it should be. 

Our representation would work also for the smoothed noise the potential —6{x) 
now being replaced by —ip * ip{x), and also in higher dimensions. But the result is 
inconsequential unless one has a detailed information on the matrix element in (4.3) for 
all N. Fortunately, the attractive (5-Bose gas can be solved through the Bethe ansatz. 
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which would still not be enough. However, for the special case {xi,...,xn) = 0, the 
formula for the propagator is sufficiently compact to proceed further. 

With this background information we return to the sharp wedge, for which 

(Z(0,f)^) = (0|e-^^*|0) (4.5) 

with shorthand |0) = |0, ...,0). To compute the propagator we use the eigenfunctions of 
H]\f, = Ex^x, which are obtained from the Bethe ansatz 

N 

^x{xi ,..., xn) = X] n , (4.6) 

P&Vn 

where the sum is over all permutations of 1,..., A^. The wave numbers Xp(a) turn out to be 
complex. They are arranged in M strings, 1 < M < N, each one containing Ua elements, 
na = N. A string has momentum eM. and Ua purely imaginary rapidities as 

A“’'’= ga + i|(n« + 1 - 2r), r = l,...,na. (4.7) 

The corresponding energy is 

M M 

riaql - ^ - ria) ■ (4.8) 

a=l ol=1 

We want to compute the generating function 

OO / inAT 

Gt{s) = ( exp [ - ^ i_Le-^(«-h/24)) t )^). (4.9) 

N=0 

The left average is for a bounded function, hence well-dehned. For the sum, however, the 
moments should not grow too fast. In actual fact 

(Z(0,t)^) (4.10) 


and the sum is ill-defined. Still, the resulting formal computation is instructive and 
nourishes the hope that one could deal with more general expectations, resp. initial 
conditions. We follow the steps of [25,53]. 

With the appropriate interpretation of the sum over A, compare with (4.13), the Bethe 
ansatz eigenfunctions are complete [54]. Hence the Wth moment can be written as 

{Z(0,()'') = 5^|(0|>l-i)|V®i‘. (4.11) 

A 


To compute the weights |(0 |\[^a)P is rather involved [55,56] and can be done systematically 
only through a proper spectral theory, see [57]. One arrives at 


Ko|%)P 


A^! det 


j(n„ -f np) + i(g« - qp) 


M 

q;,/3=1 


(4.12) 
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and the sum over eigenvalues means 


N M oo M 

n / dqa{27r)-^6{J2^c.,N) 

A M=1 ■a=lnc=l“'® o=l 


Hence 


G.fsl = 


°o (_t\M ^ r °o 

Z y-n //""p-r' Z 

M=0 ■ 0=1'^’* nc=l 


-lo«a*/24-n„ (t+s) 


X det 


\{na'+niii)+ i{qa> — qi3i) J 


M 


with the understanding that the term with M = 0 equals 1 
Using 

1 


+ n^) + i(g„ - Jq 




and a simple identity for determinants, one arrives at 

^ / _-1 \Af ^ pc>o n 

G'<(f) = E Unw n / / <' 9 P 2 ^)'‘ E 

M=0 a=ldo Jr n^=i 

Xe”“U24-nc(iig2i+s) ^ 

V Ja ',/3'=1 

Setting 7 i = (t/2)^/^, we rescale s —)■ ytSjCUa ^ diCUa, ^ da/it and obtain 


Gtilts) — det(l — PoKt sPo ), 
where the integral kernel of Kf g is given by 

« OO 

n=l 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


The projection Pq comes from the restrictions cu > 0,a;' > 0. As anticipated, the sum 
over n is badly divergent. We invoke the identity 




dyMiy^"’. 


(4.19) 


The non-uniqueness of the moment problem is reflected by the fact that Ai in (4.19) could 
be replaced by some other function out of a huge family without obstructing the identity. 
Proceeding with this choice, we set A = = 7^/2^/^ and shift y to y + + J{oj + oj') ^ 

Kt,s{oo, oo') = 2^/3 j dq{27r)-h^'^^^'-‘^^ j dyki{2^/^{y + q^ + \{uj + u')) 

OO 

X . (4.20) 

n=l 


17 






Summing over n and using the identity 


2^/3 / dg(27r)-V^^“Ai(22/3(g2 + a;)) = Ai(a; + u)Ai(a; - 


u) 


(4.21) 


one obtains 


= / dMAi(ci; + M)Ai(ci;'+ u) 


^-yt{u-s) 


X _|_ Q'ytiu-s) ’ 


(4.22) 


in agreement with Eqs. (2.6), (2.7). 

For the joint distribution at two locations, the simplicity of (4.16) is lost and one 
can proceed only under further uncontrolled approximations [58-61], which do however 
lead to the expected behavior for f —)■ oo. For flat initial conditions, instead of setting 
{xi,...,xn) = 0 in (4.3) one has to integrate over all arguments. This corresponds to 
an initial pure condensate for the (5-Bose gas which is of great interest in the context 
of quantum quenches [62]. The analogue of the generating function (2.6) is derived in 
[28]. The reader should be warned that the details are considerably more involved than 
in the example of this section. For stationary initial conditions the replica solution is 
accomplished in [26,27]. 

5 How to circumvent diverging sums 

We consider the ASEP with step initial conditions to approximate the KPZ equation with 
initial sharp wedge and introduce the modified partition function 


7 ASEP 




(6.1) 


with T = p/q, j G Z, and f > 0. As crucial observation, this family of observables satishes 
an analogue of (4.3), 


N 


((t- (5.2) 


e=i 


Here is the lattice N particle operator 


N 


^ASEP _ ^ ^ _ 


(5.3) 


£=1 


where V^/(ji,..., Jjv) = /(ji,..., j£±l,..., Jv) -/(ji, • • • Jn) with boundary conditions 
imposed through its domain, which consists of functions / satisfying 


(pV,+i + 9V/)/(J) = 0 whenever jr+i = jr + 1 


(5.4) 


In (5.2) it is assumed that ji < J 2 < • • ■ < jw- In contrast to the KPZ equation the 
“Hamiltonian” is not a Hermitian operator. 


18 



As general lore, such one-dimensional lattice systems, with only nearest neighbor inter¬ 
actions and conserved number of particles, are solvable through the Bethe ansatz [63,64], 
The reality is that most initial conditions lead to unmanageable sums over the permuta¬ 
tion group. In very special cases one has nice expressions. A particularly pleasant one is 
found in the case of step initial data, where initially all the sites {1, 2, 3,...} are occupied 
and sites {..., —2, —1, 0} are empty. As established in [65], the right hand side of (5.2) 
can be expressed through an iV-fold contour integral as 


^NiN-l)/2 

(27ri)"- 



n 

1<A<B<N 


za - Zb 
- TZb 


N 

IJexp [ 

i=l 


(Zi 


Zi{p-qy B ^ l + Zj xj, dzj 
+ l){p + qzi) ^^1 +Zi/r' Zi + T 


(5.5) 


where the integration contour C is a circle around —r, chosen with small enough radius so 
that —1 is not included, nor is the image of C under multiplication by r. It is important 
to be precise about the contours, because formulas like (5.5) really code sums of residues. 
Such contour integrals have a long history [66]. Similar formulas appear in [67] as solutions 
to (4.3). Note the important difference that unlike (4.3), the formula only holds for 
ji < j 2 < ■ ■ ■ < Jn, and there is some non-trivial work to recover moments from them. 

From the solution (5.5) one has to produce a generating function. The algebra leads 
to quantum deformations of the standard exponential function. 


eriz) = 1 /{z]t)oo 


(5.6) 


with the Pochhammer symbol {z-,t)oo = (1 ~ '2)(1 ~ zt){1 — zt'^) .... (Most common is 
the notion of g-deformat ion. It has become a r-deformation, since q is being used already 
as jump rate). [65] provide a general route to the moments and the generating functions. 
The result parallels the original derivation of exact formulas in this case by Tracy and 
Widom [42], and the formulas themselves are similar. For all complex ( except positive 
reals one has 


E(ex(Ce5<‘“*"><'“''“'®'’0.i)+)))^ = det(/ + , 


(5.7) 


where 


^ASEP 


{w, w') 


1 

27ri 



r(-i,)r(i + i,)(-c)* 


T^w — w' 


(5.8) 


with 


(5.9) 


The contour V is a. cigar around the positive integers {1,..., i?} of width 5 > 0, cutting 
the real axis between 0 and 1 and then goes straight up from R + i6 to R + ioo and likewise 
straight down from R — i6 to R — ioo, for 6 sufficiently small and i? > 0 sufficiently large so 
that for w, w' E C, a contour containing 0 and —r but not —1, and s E V, \q^w — w'\ > 0 
and \gj^t{w)/gj^t{q''w)\ < for some < oo. 

One may ask how the moment problem has been avoided. Indeed, the moments 
are expected to grow roughly like , the precise asymptotics being un¬ 
known, which is still far too fast to recover the distribution. Note that in (4.10) the 
growth is e*^^^, which demonstrates that the value of exponent for N has little to do with 
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belonging to the KPZ class. The magical point is however that because logr < 0 for 
left-finite initial data < 1. So it just comes out of this algebra. 

To complete the program, one still has to show that in the weak drive limit, as dis¬ 
cussed in Section 3.2, the generating function (5.7) converges to the generating function 
(2.6). This is a difficult piece of asymptotic analysis and still has to be accomplished. 
In the replica method such a step is avoided because the computation is directly for the 
KPZ equation with no approximation scheme. 

Acknowledgements. We are grateful for the generous hospitality at the Institute for 
Advanced Study, Princeton, where the first draft was written when both of us partic¬ 
ipated in the special year on “Non-equilibrium Dynamics and Random Matrices”. We 
thank Patrik Ferrari for help with the hgures. 


References 

[1] K. Johansson, Random matrices and determinantal processes. In: Mathematical 
Statistical Physics, Ecole d’ete Physique, Les Houches, session LXXXIII (2006), 
arXiv:math-ph/0510038. 

[2] H. Spohn, Exact solutions for KPZ-type growth processes, random matrices, and 
equilibrium shapes of crystals, Physica A 369, 71 -99 (2006). 

[3] P. L. Ferrari and H. Spohn, Random growth models. The Oxford handbook of random 
matrix theory, 782-801, Oxford Univ. Press, Oxford, 2011. 

[4] I. Corwin, The Kardar-Parisi-Zhang equation and universality class. Random Matri¬ 
ces Theory Appl. 1 (2012). 

[5] J. Quastel, Introduction to KPZ, Current Developments in Mathematics, Vol. 2011, 
125-194, Int. Press, Somerville, MA (2012). 

[6] A. Borodin and V. Gorin, Lectures on integrable probability, arXiv:1212.3351. 

[7] A. Borodin and L. Petrov, Integrable probability: from representation theory to 
Macdonald processes, Probab. Surv. 11, 1-58 (2014). 

[8] I. Corwin, Macdonald processes, quantum integrable systems and the Kardar-Parisi- 
Zhang universality class. Proceedings of the International Congress of Mathemati¬ 
cians 2014, Seoul, Corea (2014). 

[9] K. A. Takeuchi and M. Sano, Evidence for geometry-dependent universal fluctuations 
of the Kardar-Parisi-Zhang interfaces in liquid-crystal turbulence, J. Stat. Phys. 147, 
853-890 (2012). 

[10] K. Takeuchi, M. Sano, T. Sasamoto, and H. Spohn, Growing interfaces un¬ 
cover universal fluctuations behind scale invariance. Scientific Reports 1, 34, 
doi: 10.1038 / srep00034. 


20 


[11] S.G. Alves, T.J. Oliveira, and S.C. Ferreira, Universal flnctnations in radial growth 
models belonging to the KPZ universality class, Europhys. Lett. 96, 48003 (2011). 


[12] K.A. Takeuchi, Statistics of circular interface fluctuations in an off-lattice Eden 
model, J. Stat. Mech. 2012, P05007. 

[13] M. Kardar, G. Parisi, and Y.-C. Zhang, Dynamic scaling of growing interfaces, Phys. 
Rev. Lett. 56, 889-892 (1986). 

[14] L. Bertini and N. Cancrini, The stochastic heat equation: Feynman-Kac formula and 
intermittence. J. Statist. Phys. 78, 1377-1401 (1995). 

[15] L. Bertini and G. Giacomin, Stochastic Burgers and KPZ equations from particle 
systems. Comm. Math. Phys. 183, 571-607 (1997). 

[16] M. Hairer, Solving the KPZ equation. Annals Math. 178, 559-664 (2013). 

[17] M. Hairer, A theory of regularity structures. Invent. Math. 198, 269-504 (2014). 

[18] A. Kupiainen, Renormalization group and stochastic PDE’s, arXiv;1410.3094. 

[19] T. Funaki and J. Quastel, KPZ equation, its renormalization and invariant measures, 
arXiv:1407.7310. 

[20] G. Amir, 1. Corwin, and J. Quastel, Probability distribution of the free energy of the 
continuum directed random polymer in 1-|- 1 dimensions. Comm. Pure Appl. Math. 
64, 466-537 (2011). 

[21] T. Sasamoto and H. Spohn, Exact height distributions for the KPZ equation with 
narrow wedge initial condition. Nuclear Phys. B 834, 523-542 (2010). 

[22] M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. 4: Analysis 
of Operators, Academic Press, 1978. 

[23] 1. Corwin and A. Hammond, KPZ line ensemble, arXiv: 1312.2600. 

[24] A. Borodin, 1. Corwin, P.L. Ferrari, and B.Veto, Height fluctuations for the stationary 
KPZ equation, arXiv:1407.6977. 

[25] P. Calabrese, P. Le Doussal, and A. Rosso, Free-energy distribution of the directed 
polymer at high temperature. Europhys. Lett. 90, 200002 (2010). 

[26] T. Imamura and T. Sasamoto, Exact solution for the stationary KPZ equation, Phys. 
Rev. Lett. 108, 190603 (2012). 

[27] T. Imamura and T. Sasamoto, Stationary correlations for the ID KPZ equation, J. 
Stat. Phys. 150, 908-939 (2013). 

[28] P. Calabrese and P. Le Doussal, An exact solution for the KPZ equation with flat 
initial conditions, Phys. Rev. Lett. 106, 250603 (2011). 


21 


[29] P. Le Doussal and P. Calabrese, The KPZ equation with flat initial condition and 
the directed polymer with one free end, J. Stat. Mech. 2012, P06001. 

[30] J. Ortmann, J. Quastel, and D. Remenik, Exact formulas for random growth with 
half-flat initial data, arXiv: 1407.8484. 

[31] J. Ortmann, J. Quastel, and D. Remenik, A Pfaffian representation for flat ASEP, 
in preparation. 

[32] T. Sasamoto, Spatial correlations of the ID KPZ surface on a flat substrate, J. Phys. 
A, Math. Gen. 38, L549-L556 (2005). 

[33] P.L. Ferrari and H. Spohn, A determinantal formula for the GOE Tracy-Widom 
distribution, J. Phys. A, Math. Gen. 38, L557-L561 (2005). 

[34] C. A. Tracy and H. Widom, Level-spacing distributions and the Airy kernel. Comm. 
Math. Phys. 159, 151-174 (1994). 

[35] F. Bornemann, On the numerical evaluation of Fredholm determinants. Math. Comp., 
79, 871-915 (2010). 

[36] S. Prolhac and H. Spohn, The height distribution of the KPZ equation with sharp 
wedge initial condition: numerical evaluations, Phys. Rev. E 84, 011119 (2011). 

[37] P.L. Ferrari and R. Frings, Finite time corrections in KPZ growth models, J. Stat. 
Phys. 144, 1123-1150 (2011). 

[38] M. Prahofer and H. Spohn, Scale invariance of the PNG droplet and the Airy process, 
J. Stat. Phys. 108, 1071-1106 (2002). 

[39] H. Widom, On the asymptotics for the Airy process, J. Stat. Phys. 115, 1129 -1134 
(2004). 

[40] K. Johansson, Two time distribution in Brownian directed percolation, 
arXiv: 1502.00941. 

[41] 1. Corwin, P.L. Ferrari, and S. Peche, Universality of slow decorrelation in KPZ 
growth, Ann. Inst. H. Poincare B 48,134-150 (2012). 

[42] C.A. Tracy and H. Widom, Asymptotics in ASEP with step initial condition. Comm. 
Math. Phys. 290, 129-154 (2009). 

[43] P.L. Ferrari, Polynuclear growth on a flat substrate and edge scaling of GOE eigen¬ 
values, Comm. Math. Phys. 252, 77-109 (2004). 

[44] A. Borodin, P.L. Ferrari, and M. Prahofer, Fluctuations in the discrete TASEP with 
periodic initial configurations and the Airyl process, Int. Math. Res. Papers 2007, 
rpm002 (2007). 


22 


[45] A. Borodin, P.L. Ferrari, M. Prahofer, and T. Sasamoto, Fluctuation properties of 
the TASEP with periodic initial conhguration, J. Stat. Phys. 129, 1055-1080 (2007). 

[46] P.L. Ferrari, H. Spohn, and T. Weiss, Scaling limit for Brownian motions with one¬ 
sided collisions. arXiv: 1306.5095, Annals of Appl. Probab., to appear. 

[47] M. Hairer and J. Quastel, A class of growth models rescaling to KPZ, in preparation. 

[48] J. Gartner, Convergence towards Burgers’ equation and propagation of chaos for 
weakly asymmetric exclusion processes. Stochastic Processes Appl. 27, 233-260 
(1988). 

[49] H. Kesten, Aspects of hrst-passage percolation. In Ecole d’ete de Probabilites de 
Saint-Flour XIV, (P. L. Hennequin, ed.). Lecture Notes in Mathematics 1180, 
Springer-Verlag, 1986. 

[50] T. Seppalainen, Scaling for a one-dimensional directed polymer with boundary con¬ 
ditions Annals Probab. 40, 19-73 (2012). 

[51] T. Alberts, K. Khanin, and J. Quastel, The intermediate disorder regime for directed 
polymers in dimension 1-1-1, Annals of Probability 42, 1212-1256 (2014). 

[52] M. Kardar, Replica Bethe ansatz studies of two-dimensional interfaces with quenched 
random impurities, Nucl. Phys. B 290, 582-602 (1987). 

[53] V. Dotsenko, Bethe ansatz derivation of the Tracy-Widom distribution for one¬ 
dimensional directed polymers, Europhys. Lett. 90, 200003 (2010). 

[54] S.C. Oxford, The hamiltonian of the quantized non-linear Schrodinger equation, 
Ph.D. Thesis, UCLA, 1979. 

[55] P. Calabrese and J.-S. Caux, Dynamics of the attractive ID Bose gas: analytical 
treatment from integrability, J. Stat. Mech. 2007, P08032. 

[56] V. Dotsenko, Replica Bethe ansatz derivation of the Tracy-Widom distribution of the 
free energy fluctuations in one-dimensional directed polymers, J. Stat. Mech. 2010, 
P07010. 

[57] S. Prolhac and H. Spohn, The propagator of the attractive delta-Bose gas in one 
dimension, J. Math. Phys. 52, 122106 (2011). 

[58] S. Prolhac and H. Spohn, Two-point generating function of the free energy for a 
directed polymer in a random medium, J. Stat. Mech. 2011, P01031. 

[59] S. Prolhac and H. Spohn, The one-dimensional KPZ equation and the Airy process, 
J. Stat. Mech. 2011, P03020. 

[60] V. Dotsenko, Two-point free energy distribution function in (1-1-1) directed polymers, 
Journ. Physics A 46, 355001 (2013). 


23 


[61] T. Imamura, T. Sasamoto, and H. Spohn, On the equal time two-point distribution 
of the one-dimensional KPZ equation by replica, Journ. Physics A 46, 355002 (2013). 

[62] P. Calabrese and P. Le Doussal, Interaction quench in a Lieb-Liniger model and the 
KPZ equation with flat initial conditions, J. Stat. Mech. 2014, P05004. 

[63] C. A. Tracy and H. Widom, Integral formulas for the asymmetric simple exclusion 
process, Commun. Math. Phys. 279, 815-844 (2008). Erratum: Commun. Math. 
Phys. 304, 875-878 (2011). 

[64] A. Borodin, I. Corwin, L. Petrov, and T. Sasamoto, Spectral theory for interacting 
particle systems solvable by coordinate Bethe ansatz, arXiv:1407.8534. 

[65] A. Borodin, I. Corwin, and T. Sasamoto, From duality to determinants for q-TASEP 
and ASEP, Annals Probab. 42, 2314-2382 (2014). 

[66] G.J. Heckman and E.M. Opdam, Yang’s system of particles and Hecke algebras. 
Annals Math. 145, 139-173 (1997). 

[67] A. Borodin and I. Corwin, Macdonald processes, Probab. Theory Relat. Fields 158, 
225-400 (2014). 


24 


